A NONLINEAR MODEL DESCRIBING A SHORT WAVE LONG 
WAVE INTERACTION IN A VISCOELASTIC MEDIUM 
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Abstract. In this paper we introduce a system coupling a nonlinear Schrodin- 
ger equation with a system of viscoelasticity, modeling the interaction between 
short and long waves, acting for instance on media like plasmas or polymers. 
We prove the existence and uniqueness of local (in time) strong solutions and 
the existence of global weak solutions for the corresponding Cauchy problem. 
In particular we extend previous results in [Nohel et. al., Commun. Part. 
Diff. Eq., 13 (1988)] for the quasilinear system of viscoelasticity. We finish 
with some numerical computations to illustrate our results. 



1. Introduction 

In \T , D.J. Benney initiated the study of mathematical models describing the 
interaction between short and long waves in fluids, namely capillary and gravity 
waves or internal and surface waves. This has been developed by many authors, 
such as M. Tsutsumi and S. Hatano in the pioneering papers [211 [25], and more 
recently, in the framework of quasilinear systems, in [TUl [HI [121 [SI E] and (15] , 
Their numerical study was initiated in [I] and developed in |2]. 

In this paper, we will apply the ideas of these authors to the case of viscoelastic 
fluids, arising, for instance, in plasma physics (in the study of helioseismology, cf. 
[6]), or in the study of polymers |18j . To this purpose, we introduce a model coupling 
the nonlinear Schrodinger equation (modeling the short waves) with a quasilinear 
system of viscoelasticity (modeling the long waves) . This last system was studied 
in [20]. Thus, we analyze the following Cauchy problem. 



(1.1) 



iut + Uxx = uv ^ a\u^u 

wt = i<jiv))x + i\u\^)x + f k{t~T)[{aiv))x + {\u\%]{;T)dT, 

JO 



where a; € M, t > 0, a is a real constant, and i = ^/—l. Here, u{x,t) e C is the 
envelope of the short waves, v{x, i) € M is the deformation gradient, ■w{x, t) £ R the 
velocity of the long waves, a G C'^(M) is the stress function verifying a' > (Tq > 



(hyperbolicity), and k is a given kernel. We consider (1.1 1 with the initial data 
(1.2) Uo,vo,WoeH\R). 
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The integral term in the right-hand side of ( |1.1[ ) represents the memory effects due 
to the viscoelastic structure of the fluid. 

We now recall a transformation due to R.C. MacCamy (see also PD|): Let 
q{t) be the resolvent kernel associated with k, i.e., q is the solution of the linear 
Volterra equation 



(1.3) 



q{t)+ / k{t - T)q{T) dT, 



r > 0. 



Convolving the third equation in (1.1) with q{t), it is not difhcult to obtain for 
smooth solutions 



(1.4) 



k{t-T)[{a{v))^ + {\u\'^)^]{x,T)dT ^ / q{t-T)wtix,T)dT 



q{0)w{x,t) — q{t)wo{x) + / q' {t — t)w{x,t) dr 



Thus, for smooth solutions, (|l.l|)(|1.2|) is equivalent to the Cauchy problem 
(1.5) 



< 



lut + Uxx = uv + a\u\ u 
vt = Wx 

[wt = {cr{v))x + {\u\^)x + Hw), 



X e 



t > 0, 



u{x,{)) — uq{x), v{x,0) — vq{x), w{x,Q) — wq(x), 



3'{w){x,t) — q{0)w{x,t) — q{t)wo{x) + / q' {t ~ t)w{x,t) dr. 



(1.6) 
where 

(1.7) 

In the case where u = 0, that is, the nonlinear viscoelasticity system, the ex- 
istence of a global (in time) weak solution for initial data in L°° (R) n (M) and 
k € C^([0, -l-oo)) was proved in [5D] (see also [S] for a different model), by the van- 
ishing viscosity method applied to both variables v and w and the compensated 
compactness method [531 [TB] based on ideas introduced by C. Dafermos in [7] (see 
also [TT], and for related results p.[TU]). 

In this paper, we start by proving in Section [2] the existence and uniqueness 
for local (in time) strong solutions for ( 1.5 ),( |L6| ) with initial data {uq,vo,wo) € 
H^(M.) X (ff2(K))^ by applying a variant of Kato's theorem 07]. In Section [s] and 
for a special class of stress functions a, we apply the physical vanishing viscosity 
method (that is, only in the velocity variable w) and a variant of the compensated 
compactness method introduced by D. Serre and J. Shearer in [35] (see also [T2^) to 
obtain the existence of a global (in time) weak solution for (1.5), (1.6) with initial 
data {uo,vo,wo) G (i7^(M))'^. Finally, in Section |4] we present some numerical 
simulations to illustrate the behavior of the solutions in a special case. 

2. Local (in time) existence of strong solutions 
Let uo e H^{E.), vo € H'^{M.), wq E H^(. 



(1.5),(1.6), we introduce the Riemann invariants 

l = w+ \/<^'iOdL r = w- 
Jo 

We derive I — r = f{v), for some / one-to-one and smooth, and w 



To study the Cauchy problem 



_ l+r 
2 ■ 
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For classical solutions the Cauchy problem (1.5), (1.6) is equivalent to 



(2.1) 



iut + — uv + a\u\ u 
h-V^l. = {\u\% + ^J{l + r) 

rt + V^r, = {\u\^)x + + r), 



by setting v = f (/ — r) = v{l,r) and with 3^ given by (1.7). We take as initial 
data m(-,0) = uq e H^{R), l{-,0) = Iq e H'^{R), r(-,0) G H'^{R), with 



(2.2) 



Wq 



To obtain a local strong solution of the Cauchy problem (2.1), (2.2), we consider, 
using the technique employed in |21j and |14| . an auxiliary system with non-local 
source terms. This is necessary in order to write the system (2.11,(2.2) without 
derivative loss (see [H]). Thus, we consider the following system, 

1 



(2.3) 



iFt + 



Fv + au^F + 2a\urF - 



zu{lx + r^) 



lt~V^lx = {\u\% + \ni + r) 



1 



where F is the complex conjugate of F and F, u are defined by 

u(x, t) — Uq{x) + 



/ f (x, s) ds, 
Jq 



u{x, t) = (A - iy^{a\u\'^u + u{v - 1) - iF) 
(see [13] for the motivation behind this definition). The initial data are 
(2.4) F{-,0)^FoeH\R), l{-,0) ^ lo e H^iR), 



with Iq, vq given by (2.2). 

We will prove the following result: 

Theorem 2.1. Let {Fq, Iq, tq) e x x . Then, there exists T* > (depend- 
ing on {FQ,lQ,ro)) such that for all T < T* there exists a unique solution (F,l,r) 
of the Cauchy problem (2. 3), (2. 4) with 



{F,l,r) e C^i[0,T];H'-'^) x {[0,T]; H^'^) x i[0,T]; H^-^), 



J =0,1. 



From this result and from the definitions of F, u, reasoning as in ^A, Lemma 
2.1] (see also [H]), it is easy to derive the following result for the system ( |1.5| ,( 1.6 ): 

Theorem 2.2. Let {uq,vq,Wq) e x x H'^ . Then, there exists T* > 
(depending on {uq,vq,wq)) such that for all T < T* there exists a unique solution 
{u,v,w) of the Cauchy problem (1.5), (1.6) with 

{u,v,w) e C^([0,T];i/3-2i) x {[0,T]; H^'^) x {[0,T]; H^-^), j = 0,1. 
Proof of Theorem\2.1\ In order to apply a variant of Kato's well posedness result, 



Theorem 6 in [T7], we put the Cauchy problem (2.3),(2.4| in the framework of real 
spaces by introducing the new variables Fi — diF, F2 — 3F, ui = diu, U2 — 3w. 
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By s e ttin g U = {Fi,F2,l,r) and Fio = ^Fq, F20 — 3i^o, the Cauchy problem 
(2. 3 1,(2.4) can be written as follows, 



(2.5) 
where 



A{U) 



( Ut + AiU)U = g{t,U) 
jc/(-,0) = C/o, 

A 

-A 

-^cr'(u)9^ 

^a^v)d. 



9it,U) = 



2a\u\^F2 - a{ul - u^)F2 + 2auiU2Fi + F2V + \u2{lx + r,^) 
2a\u\'^Fi - a{u\ - ul)Fi - 2auiU2F2 - Fiv - \u2{lx + rx) 
(|u|2), + lj(; + r) 

(|un, + ij(; + r) 

and C/o = (i^io, ^20, '0, ^-q) e y = {H^{M.)f x {H'^{R)f. Note that the source term 
is non-local. 

In what follows we use the notations of Tf, paragraph 7]. Set X — (_ff^^(M))^ x 
(L^(M))^ and S = {1 — A) Id, which is an isomorphism S :Y X. Furthermore, 
we denote by Wr the open ball in Y of radius R centered at the origin. 

We need to check several assumptions in order to apply Theorem 6 of |17j. First, 
it is necessary that the semigroup generated by the operator A above verify 

(2.6) ||e-*^(^)|| < e"*, 

for some real oj, for alH > and U £ Wr. Observe that it is enough to prove this 
only for the operator 

' 



a{l,r) 







since the remaining part of A corresponds to the Schrodinger (contraction) group. 
But (2.6 1 is proved in [171 paragraph 12], with u given by 

a; = 1 sup \\d^a{l,r)\\ < c{R), 

with c : [0, +00) — !■ [0, +00) continuous. 

Next, we must check that for U e Wr, the property SA(\J)S^^ = A{U)+B{U) 
is valid for some B G -C(X). This is proved in 17, paragraph 12]: for {I, r) in a ball 
W of (i/2(M))2^ we have 

(1 - A)a{l, r)(l - A)-i = a{l, r) + So(/, r). 



where 
and 



Bo(;,r) = [(1 - A),a(?,r)](l - A)-i e L{{L^mf ) 



B{U)^ 



with 



Bo(Z,r) 

denoting the matrix commutator operator. 
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Now, consider a pair £ C{[0,T];Wr), U = (Fi, F2, Z, r), [/* = {F^ , ,1* ,r*[ 

It is easy to see that g verifies, for fixed T > 0, \\g{t, U{t)\\Y < C{R,T), t e [0,T], 
iiU e C{[0,T];Wr). We obtain, in the same way as [II ITS], 

(2 7) \\9{-,U)-g{-,U*)\\LHo.T',x)<c{T') sup \\U{t) - U* {t)\\x 

^ ' 0<t<T' 



where < T' < T and c(T') is a non-decreasing continuous function such that 
c(0) = 0. Finally, applying Theorem 6 in [17], replacing the local condition [ITl 
equation (7.7)] by (2.7), we obtain the result. This completes the proof of Theo- 
rem [HI □ 

3. Global existence of weak solutions for a class of stress 

functions 

Now we will consider the question of global (in time) weak solutions of the 



Cauchy problem (1.5),(1.6|. For a special class of stress functions a, we will obtain 
an extension of the result in [5D] for the system of nonlinear viscoelasticity and 
in |12j for the system of nonlinear elasticity coupled with the nonlinear Schrodin- 
ger equation. We employ the adaptation of the compensated compactness method 
developed by D. Serre and J. Shearer in 22 for the system of nonlinear elasticity, 
which extends earlier results of L. Tartar i23j and R.J. DiPerna P^. 

Following we define S(w) — a{v) d£, (we may assume cr(0) = 0) and we 
impose the following conditions on the stress function a G C'^(M): 

HI: <j'{v) > (JQ > for some constant (Tq. 

H2: cr"(Ao) = and ct"(A) ^ for A ^ Aq. 

rr" rr'" rr" rr'" 

H3: ^^ei2(M), — eL°°(M). 



H4: 



S(«) 



as - 



00, and there exist constants c > and g > 1/2 such that 



In particular, we have S(f) > 



To simplify, we also suppose that v £ 



implies '^(v) dx < 00. A typical example is given by a{v) = + v. 

We now follow the ideas developed in [T^ and introduce a physical viscosity 



approximation of the Cauchy problem (1.5), (1.6). For small e > 0, we consider the 
system 



ilUt 
vt- 
Wt 

with initial data 
(3.2) 

Observe that although in 



uv + a\u\ u 



X e 



t > 0, 



{a{v))^ + {\u\'^),^ + Jiw) + ew^ 



uq, vq, Wo e ff^(M). 



in L° 



the authors consider initial data for which vq, wq are 
only, this is not possible in our case, since the method of Serre and Shearer 



requires initial data with stronger regularity. Of course, this does not exclude the 
formation of discontinuities for t > 0. 

The proof of the following lemma is an easy adaptation of the first part of the 
proof of Lemma 6.1 in [T^], and so we omit it. 
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Lemma 3.1. Let [u, v, w) G (C([0, T); i7^(M))^, T > Q, be a solution to the Cauchy 
problem (3.11,(3.2) for fixed e > 0. We have in [0,T] 



(3.3) 



(3.4) 



^ / \u\^dx = 0, 



d 
dt 



[ I \ux\^dx+^ I \u\'^dx+ I v\u\^dx+\- f w'^ dx + f J:{v) dx] 
Jr ^ Jm Js. 2 



3^(10)10 dx + e w'^dx — 0. 
; Jm 

Now, we recall from 2D] the inequality 



(3.5) 



J{w)wdxdT < C{T)(l + 



' dx dr 



where C(-) is a continuous function in [0, +00), independent of e. 

We now prove some estimates which will be needed to prove the well-posedness 
of the viscous Cauchy problem (3.11,(3.2). 



By the Gagliardo-Nirenberg inequality we have, by (3.3), 

|u|'*da; < Co||w||i2||u2,||i2 = Co||uo|li2 



and since E{v) > we derive from (3.4), (3. 5) for < i < T, 



(3.6) 



w'^{x,t)dx < C{T) + C{T) 



w'^{x, t) dx dr 



and so, by the Gronwall inequality and (|3.4|)-(|3.6|) we get for < i < T, 
(3.7) 



|Ma;(a;,i)| dx + \w{x,t)\ dx + 'S{v{x,t))dx 



if [ wldxdr < C(t), 
Jo Jm 



(3.8) 



J{w)wdxdT < C{T), 



for some continuous function C on [0, +00) independent of e. 
Next, we will prove the following estimate, for e < 1, i e [0,T]: 
ft f 

a'{v){v^fdxdT + e^ I {v^fdx<C{T), 



(3.9) 

where C is a continuous function on [0, +00) independent of e. For this purpose 
we follow the ideas in [22^ proof of (8)] (see also [H]). We deduce from (3.1), for v 
smooth enough. 



WtVx - a' (v){vx)'^ dx ^ / {\u\'^)3:Vxdx + / 3'{w)vxdx + t / w^xVxdx 



and 



d 
"dt 



Wxvdx+ / {wx)'^ dx — / a' {v){vx)'^ dx 



\u\'^)xVxdx + / 3'{w)vx + 



e d 
2dt 



(vx)"^ dx 
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(since, for w smooth enough, — ^ J^WxV dx — — /g Wxtv dx—J^^ WxVt dx = WtVx dx- 
/jjj WxVt dx and vt — Wx)- 

Integrating over (0,t) we obtain, with ^0(2;) = v{x,0), wo{x) ~ 'w{x,0), 



Wxvdx+ / woxVodx+ / / {wx)'^dxdT 



{v){vxY dx dr 



i\uf)xVx dxdr - 



From — J^WxV dx = wvx we derive 



3^{w)vx dxdT+ ^ I {vx f - {vox f dx. 



C I a\v){vxf 




dxdT+- (vx)^ dx < - j {vxY dx -\- 



1 



dx 



(3.10) 



|it;oa;||wo| dx 



{vqxT dx + e I \{wxYdxdT 



\u\\ux\\vx\ dx dr + 



3^{w)vx dx dr 



Moreover, we easily get from (3.71, for a fixed (5 > 0, 



(3.11) 



dxdr < C{5)C{T) + 



<^o Jo 



a'{v){vx)'^ dx dr, 



and also by the definition of 3^, 



(3.12) 



3^{w)vx dx dr 



< S 



{vxfdxdT + C{6)C{T) / / w^dxdT 



< 



f^o Jo 



cj'{v){vxfdxdT + C{5)C{T) 



' dx dr. 



From (3.10 )-(3.12 1 we derive, choosing 5 = ^ and multiplying by e, the estimate 
dSj] ) for t G (o,r]. 

Let us now analyze the problem of the existence of a unique solution (u, w, w) G 
(C([0, +00); 7Ji))3 of the Cauchy problem ([3l]) , ([3|2| for fixed e G (0, 1]. We assume 
e = 1 without loss of generality. Let us begin with the study of the existence and 
uniqueness of a local (in time) solution. Let < T < +co and introduce the complex 
Banach space (not to be confused with the spaces introduced in Section |2]) Xt ~ 
C([0,T];iJi) and the real Banach space Yt = C([0, T]; i/^), with the standard 
norms. Let us consider the product space x {B^Y where = {tt G Xt ■ 
\\u\\xr < R} and (5^)2 = {{v,w) G : \\v\\y^ < R, \\w\\y^ < R}. 

Given {u,v,w) G x [B]^)'^ we consider the mapping 

{u,v,w) I— ^{u,v,w) = {u,v,w) G Xt x (1t)^ 
where (w, v, w) is the unique solution of the linear problem 

iut + Uxx — a\u\^u + vu, u(0) = uq, 

Wt - Wxx = {(^{v))x + {\u\'^)x + 3'(w), ^(0) = Wq, 

v{t) = vo+ Wx dr, 



(3.13) 
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with uo,vo,wo e and such that w € L^{0,T; H^),wt e L'^{0,T] L'^). With the 
usual method of semigroups, we have 

u{t) = e'^'uo - i f e^^'-'\a\u\'u + iu){s) ds 



and 

(a(t))), + (|un,+J(^)eC([0,T];L2) 
with the estimate [20l (3.42)] 



(3.14) / / \3:{wi)-J{w2)? dxdT <C{T) I I \wi - W2\^ dx dr. 



Reasoning as in [20', proof of Theorem 3.1], we can, for convenient T and R > 
maxdwolffi, |wo|hi, IwqIhOi ^PPly the Banach fixed point theorem for strict con- 
tractions to obtain existence and uniqueness of a local (in time) solution of the 
viscous system (3.1), (3. 2). Moreover, we have 

Wt - W^.^ = {a{v))r, + (|up):r + J(w) 



and, from the previous estimates (3. 7), (3.91, (3.14) (with lii = w, 11)2 = 0) and 
(3.3), we derive the a priori estimate 

\\Wt - W^^\\l2(^o^T;L2) < c{T), 

with c G C([0, +oo[; M_|_), and this implies, by standard regularity estimates for 
parabolic equations, that w G L^{Q,T; H^) and |j w|ji2(Q ^^.^2-) < c(T). This in turn 
gives similar a priori estimates for ||wt||x,2(o, t^l^) and HwUcQo.Tliffi)- 

We are now in a position to pass to the limit e — > and state the main result in 
this section. 

Theorem 3.2. Assume {uo,vo,wo) e (iJi(M))3 and let Lf^^{R x [0,+cx))) be the 
space of functions v measurable in M x [0, +00) such that 

I](w) dx dt < 00 for each compact K dR x [0, +00). 

IK 

Then there exists 

{u,v,w) e ij^,((0,+(X));i?i) X Ll^(R x [0,+oo)) x LI^{R x [0,+^)) 
such that 

00 



u6t dx dt — / / UxOj.dxdt + I uq{x)9{x,Q) dx 
Jo Jr Jr 

vu9 dx dt + a / / \u\^u9dxdt, 



10 JR Jo Jr 

for all complex-valued € Cq (M x [0, +00)), and 



{vcfit — w(j)x) dx dt + j vo{x)(/){x,0) dx 
Jr 

wipt — criv)4'x + 3^{w)ip dx dt 
Wq{x)iP{x,0) dx -\- / / {\u\'^)x^jJ dx dt — 0, 



for every real-valued (f),ip & Cq{M. x [0, +00)). 



SHORT WAVE LONG WAVE INTERACTIONS IN A VISCOELASTIC MEDIUM 



9 



Proof. We follow the ideas in [22l[T2] (see also [4l): for each e e (0, 1] let {ue, Ve, We) G 
(C([0,+oo);i7i)) be the unique solution of the Cauchy problem (|3.1|,(|3.2|). From 
the system (3.1), (3.3), and (3.71, we derive, for fixed T > 0, 

{uje bounded in i°°((0, T); iJ^), 

{uei}e bounded in L°°((0, T); i?'^), 

and so, by a well known compactness result, {u^} is in a compact set of ^^(o, T; L'^{Ir)) 
for each interval Iji = {—R,R), R > 0. Furthermore, there exists u S H^{M.) such 
that (for a subsequence) ^ u in ij^^(0, +oo; iJ^(M)) when e — > 0. By a stan- 
dard diagonalization method, we conclude that in fact, for a suitable subsequence. 

We also have, by (3.7), 

{we}e bounded in ^^^^(M x [0, +oo)), 
{fe}e bounded in L^(M X [0,+oo)). 



By ( |3.7[ ) and ( |3.9[ ) we derive 
ft 



(3.15) 



{w,^Y + a\v,){v,^Ydxdt<C{T) 



where G(t) is a continuous function depending the norm of the initial data, but 
not on e. 

Now consider the quasilinear system 



(3.16) 



Vt = 



and, for (u, w) G M^, let ij^v, w), q{v, w) be a pair of smooth convex entropy-entropy 
flux pair for (3.16) such that rj^, rj^^ and ri^wl\f^ are bounded functions in M^. 
From (3. 3), (3. 7), (3. 15) and since, for each interval = [— i?, i?], i? > 0, we have 

J{w^)dxdt < Cfl(t), 

JIr 

where 0^(1) is a continuous function on [0,-|-oo) independent of e, we derive from 
the system (3.1) (cf. similar estimates in [22111] and [12]) that 



dtri{Ve, We) + d^qiVe, We) 



belongs to a compact set of W, 



-1,2/ 



loc 



X [0, -|-cxd)). Hence, in view of the assumptions 

x[0,+cx))))' 



H1-H4, we may apply the result on compensated compactness of D. Serre and 
J. Shearer [22] to conclude that {(ve,We)}e is pre-compact in {LI^J^ ' ''^ 
Hence, there exist a subsequence {(ue,Ve,We)}e and 

{u, V, w) e LJ^,(0, +oo; H^) x Lf^^{M. x [0, +oo)) x x [0, -|-c3o)) 

such that 

{ue, Ve,We) — )■ (u, V, w) in (Lj'o(,(IR X [0, +00)))^, e 0. 

Now take functions 9, (j) and tjj as introduced in the statement of Theorem |3.2[ 
multiply the first equation of system (3.1 1 by 0, the second one by (j) and the third 
one by tfj and integrate over M x [0, +oo). After integrating by parts, we pass to 
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the limit as e — > in view of the above convergences. This allows us to prove that 
(u, ?j, w) are weak solutions of the system ( 1.1 ), since we have for some T = T{ip), 

Wexx'ipdxdt <c^e^^^(e / / (wex)'^ dx dt] 

^ .In ./tu 



1/2 



which goes to zero as e — > by (3.7). This completes the proof of Theorem 3.2 □ 



4. Numerical experiments 



In this section we present some numerical experiments on the system (1.1), to 



illustrate our results. We use the more convenient formulation ( 1.5 )-( 1.7). Accord- 
ing to [5D], if the kernel k in ( 1.1 ) is the derivative of a C^, positive, decreasing and 



convex function on [0,-|-oo), then the term q{0)w{x,t) in J has a damping effect 
and the solutions should remain classical and globally defined. 

Thus, we will consider the particular case k{t) = e^* in (1.1), which gives, after 
some elementary calculations, q{t) = e~^* in ( |1.5| ). 

The stress function is taken as (j{v) = v^ + v, which verifies the conditions H1-H4 
111 Section [3] The memory term in (|1.5[) reads 



(4.1) 



3^{w) = w{x, t) 



^wq{x) 



e2(''-*)u;(x,s) ds. 



4.1. The numerical scheme. We use a fourth-order explicit Runge-Kutta scheme 
for the time-stepping along with standard finite difference discretizations of the 
space derivatives. As is usual in the simulation of problems posed on the whole 
line, we restrict ourselves to a bounded domain [Li, L2] and to initial data decaying 
exponentially for large |a;|, and set the boundary conditions to zero. 

Thus, we are given a spatial mesh size h > 0, J = (L2 — Li)/h (which we suppose 
is a natural number without loss of generality), a time step r, a suitable discretiza- 
tion of the initial data, (uoj, woj , w'oj)j=o,...,Jj an approximation (w", w", w")j=o,...,j 
of {u{xj,tn), v{xj,tn), w{xj,tn)) at time i„ = nr, n G N, and Xj = Li +jh. We then 



obtain (u^ 
equations 



n+l „,n+l 



,W 



by solving by the fourth-order Runge-Kutta method the 



(4.2) 



dt ^ 



dt ^ 



(^^;Vl-K + "^l)-™"KI' 



«+l 



-1) 



1 

' 2h 

= ^(-K-Vi)--K--i)) + iK+i 



2e 



where 5^, the discretization of the non-local integral term in (4.1 ) is defined recur- 
sively by 



Note the introduction of a viscosity term in the third equation on (4.2 ) to improve 



stability, as is customary in the numerical treatment of hyperbolic equations. 
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Figure 1. The initial data (4.3) 



As initial data (see Figure [ij, we have set 

uoix) = Ce^°"cosh~\V50x), 

(4.3) uo(a;) = Ccosh~^(\/20(a;-0.1)), 

woix) = Ccosh"i(V20(x + 0.1)). 

For other simulations of short wave long wave interactions in the case of a Schro- 
dinger equation coupled with a nonlinear conservation law, see ^[l]. 

4.2. Numerical results. We now present some results obtained by this numerical 
method. In Figure[I] we present the initial data ( |4.3[ ). Next, In Figures|2[|5j we plot 
the numerical solution, respectively, at times T = 0.001, T = 0.01, and T = 0.1, 
with a close-up of the interaction region in Figure [5] 

We can clearly observe in these simulations the interaction between the short and 
the long waves. In particular, there is the formation of a train of small frequency 
waves in the v and (more strongly) in the w variables which can clearly be discerned. 
These new waves, due to the nonlinear coupling, are coherent with the oscillations 
in the Schrodinger variable u, as can be seen in Figures [3[|5j 

Thus, our simulations show creation of new waves due to the interaction effects 
and thus allow us to see a new qualitative property of the solutions to the system 
under consideration. 
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dt = 5*10"-6, 5500 points, dx= 0.0055, T = 0.001 
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Figure 2. Computed solution of system (1.5 1,(1.6 1, T = 0.001. 



dt = 5*10"-6, 5500 points, dx= 0.0055, T = 0.01 




Figure 3. Computed solution of system (1.5), (1.6), T = 0.01. 
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